A.

      B.

Figure 1. Stem density distributions by cycle in Richmond National Battlefield Park (RICH). Cycle 1: 2007 – 2010; Cycle 2: 2011 – 2014; Cycle 3: 2015 – 2018; Cycle 4: 2019 – 2022; Cycle 5: 2023. Note that cycle 5 is only partially complete. Figure A shows regeneration densities by size class. Figure B shows tree density by diameter at breast height (DBH) increments. Error bars are 95% bootstrapped confidence intervals to account for non-normal error. An asterisk denotes a linear DBH distribution that is indicative of long-term recruitment failure.



Figure 2. Regeneration Debt status for the most recent 4-year census (2021 – 2023) in Richmond National Battlefield Park (RICH). Flat Tree Diam. Dist. stands for Flat Tree Diameter Distribution. If TRUE, the density of small trees is lower than expected due to chronic regeneration failure. Sapling/Seedling Composition is the % of total stems composed of native canopy-forming species. Sorenson is a measure of how closely the regen. layer matches canopy composition. For more details on metric calculations, thresholds and assessment of status, see Miller et al. 2023.


Figure 3. Deer browse impacts by cycle in Richmond National Battlefield Park (RICH). Cycle 1: 2007 – 2010; Cycle 2: 2011 – 2014; Cycle 3: 2015 – 2018; Cycle 4: 2019 – 2022; Cycle 5: 2023. Note that cycle 5 is only partially complete.


Figure 4. Loess smoothed changes in stem density for seedlings (A) and saplings (B) by species group and sample year in Richmond National Battlefield Park. Specified loess span approximated a linear trend between consecutive sample events per panel.



Figure 5. Loess smoothed changes in tree stem density (A) and basal area (B) by species group and sample year in Richmond National Battlefield Park. Specified loess span approximated a linear trend between consecutive sample events per panel.



Figure 6. Loess smoothed changes in invasive plant percent cover by guild and sample year in Richmond National Battlefield Park. Specified loess span approximated a linear trend between consecutive sample events per panel.



Table 1. Average plot-level seedling and sapling stem densities (stems per m2) and stocking index (at the 2m radius scale) by cycle. Only native, canopy-forming species are included. Note that Fraxinus spp. (ash species) are no longer considered canopy-forming species. Cycle 1: 2007 – 2010; Cycle 2: 2011 – 2014; Cycle 3: 2015 – 2018; Cycle 4: 2019 – 2022; Cycle 5: 2023. Plots are sampled in 4-year rotations, with a quarter of the plots sampled each year (i.e. a panel). The stocking index quantifies whether current regeneration densities are sufficient to restock a forest canopy. The index is a weighted sum of seedling and sapling densities where larger seedling size classes get higher weights. Cells highlighted in green are plots that meet the minimum management target of 0.25 seedlings/m2 and 0.14 saplings/m2 or have a stocking index >100. These are the same thresholds used in Figure 2 and are described in Miller et al. 2023. Note that cycle 5 is only partially complete.

table output
Seedlings per sq.m
Saplings per sq.m
Stocking Index
Plot Panel Cycle 1 Cycle 2 Cycle 3 Cycle 4 Cycle 5 Cycle 1 Cycle 2 Cycle 3 Cycle 4 Cycle 5 Cycle 1 Cycle 2 Cycle 3 Cycle 4 Cycle 5
15 1 0.00 0.42 0.25 0.08 0.08 0.07 0.06 0.02 0.02 0.02 14.83 21.11 10.56 8.42 8.42
16 1 0.00 0.00 1.33 0.25 2.25 0.02 0.00 0.01 0.02 0.07 0.00 0.00 17.84 10.56 67.23
19 1 0.08 0.17 0.00 0.00 0.08 0.00 0.00 0.00 0.00 0.01 2.14 2.14 0.00 0.00 8.42
20 1 0.58 0.25 0.25 0.25 0.08 0.04 0.01 0.02 0.02 0.02 22.12 3.14 10.56 11.56 9.55
63 1 0.50 0.00 0.17 0.00 0.08 0.01 0.01 0.01 0.01 0.00 10.43 0.00 2.14 0.00 1.01
64 1 0.08 0.00 0.00 0.00 0.08 0.01 0.00 0.00 0.00 0.00 1.01 0.00 0.00 0.00 1.01
72 1 0.33 0.00 0.00 0.17 0.25 0.33 0.07 0.06 0.01 0.00 192.89 7.41 7.41 2.14 3.14
73 1 0.42 0.00 2.10 0.75 0.50 0.21 0.17 0.12 0.07 0.06 27.52 29.66 36.32 11.56 7.29
129 2 0.25 0.33 0.17 0.33 - 0.09 0.13 0.08 0.08 - 34.81 35.94 17.97 20.11 -
131 2 0.17 0.08 0.25 0.67 - 0.00 0.00 0.00 0.00 - 2.14 1.01 3.14 8.42 -
133 2 0.33 0.25 0.25 0.17 - 0.32 0.25 0.26 0.14 - 81.43 69.62 62.20 11.56 -
134 2 0.00 0.00 0.00 0.00 - 0.01 0.01 0.01 0.01 - 0.00 0.00 7.41 0.00 -
148 2 0.25 0.25 0.00 0.08 - 0.13 0.07 0.05 0.02 - 85.07 17.97 7.41 1.01 -
153 2 0.25 0.08 3.42 2.00 - 0.04 0.04 0.18 0.27 - 23.00 8.42 148.53 306.75 -
154 2 0.33 0.58 0.58 1.92 - 0.25 0.22 0.21 0.19 - 86.71 76.15 61.32 76.91 -
155 2 0.00 0.00 0.00 0.00 - 0.00 0.00 0.00 0.00 - 0.00 0.00 0.00 0.00 -
203 3 0.33 0.33 0.17 0.08 - 0.18 0.25 0.27 0.28 - 254.47 298.83 143.88 193.15 -
207 3 0.50 0.17 0.08 0.08 - 0.02 0.05 0.04 0.02 - 8.42 16.96 9.55 9.55 -
215 3 0.08 3.17 2.08 0.58 - 0.15 0.19 0.20 0.15 - 24.38 96.76 145.27 60.19 -
219 3 0.67 0.67 1.17 - - 0.11 0.14 0.14 - - 57.93 51.65 43.23 - -
226 3 0.08 0.17 0.33 0.17 - 0.01 0.01 0.04 0.04 - 8.42 11.56 28.53 25.38 -
230 3 0.50 0.75 0.58 0.42 - 0.05 0.06 0.05 0.06 - 29.28 42.98 18.98 25.26 -
233 3 0.75 0.33 1.00 0.58 - 0.06 0.07 0.06 0.05 - 46.24 77.79 79.80 22.12 -
238 3 0.00 0.00 0.00 0.00 - 0.05 0.05 0.01 0.01 - 22.24 22.24 0.00 0.00 -
303 4 0.25 0.42 0.42 0.58 - 0.06 0.06 0.06 0.06 - 17.97 20.11 13.70 14.70 -
304 4 0.42 0.92 0.92 2.00 - 0.14 0.13 0.11 0.09 - 15.83 55.54 35.56 49.26 -
306 4 0.00 0.00 0.00 0.00 - 0.07 0.07 0.06 0.01 - 14.83 7.41 7.41 0.00 -
307 4 0.25 0.42 0.17 0.17 - 0.02 0.01 0.05 0.06 - 5.28 6.28 25.38 32.80 -
309 4 0.00 0.08 0.00 0.00 - 0.02 0.01 0.00 0.00 - 14.83 2.14 0.00 0.00 -
311 4 0.17 0.17 0.17 0.17 - 0.02 0.02 0.01 0.01 - 2.14 2.14 2.14 2.14 -


Table 1 (cont.). Average plot-level seedling and sapling stem densities (stems per m2) and stocking index (at the 2m radius scale) by cycle. Only native, canopy-forming species are included. Note that Fraxinus spp. (ash species) are no longer considered canopy-forming species. Cycle 1: 2007 – 2010; Cycle 2: 2011 – 2014; Cycle 3: 2015 – 2018; Cycle 4: 2019 – 2022; Cycle 5: 2023. Plots are sampled in 4-year rotations, with a quarter of the plots sampled each year (i.e. a panel). The stocking index quantifies whether current regeneration densities are sufficient to restock a forest canopy. The index is a weighted sum of seedling and sapling densities where larger seedling size classes get higher weights. Cells highlighted in green are plots that meet the minimum management target of 0.25 seedlings/m2 and 0.14 saplings/m2 or have a stocking index >100. These are the same thresholds used in Figure 2 and are described in Miller et al. 2023. Note that cycle 5 is only partially complete.

table output
Seedlings per sq.m
Saplings per sq.m
Stocking Index
Plot Panel Cycle 1 Cycle 2 Cycle 3 Cycle 4 Cycle 5 Cycle 1 Cycle 2 Cycle 3 Cycle 4 Cycle 5 Cycle 1 Cycle 2 Cycle 3 Cycle 4 Cycle 5
314 4 0 0 0 0 - 0.48 0.41 0.28 0.22 - 111.09 81.43 51.9 44.48 -
317 4 0 0 0 0 - 0.02 0.02 0.02 0.01 - 0.00 0.00 0.0 0.00 -


Table 2. Average % cover and number of invasive species detected in each plot by cycle. Cycle 1: 2007 – 2010; Cycle 2: 2011 – 2014; Cycle 3: 2015 – 2018; Cycle 4: 2019 – 2022; Cycle 5: 2023. Plots are sampled in 4-year rotations, with a quarter of the plots sampled each year (i.e. a panel). Note that cycle 5 is only partially complete. Plots with at least 10% average cover or 4 or more invasive species present are in yellow. Average cover and species counts may be increasing over time due to protocol changes, including additions to the indicator species list over time, a plot search for indicator species starting in 2009 and all woody species being sampled in quadrats starting in 2019. More details can be found in the MIDN Summary of Major Protocol Changes

table output
Average % Cover
Number of Species per Plot
PlotCode PanelCode Cycle 1 Cycle 2 Cycle 3 Cycle 4 Cycle 5 Cycle 1 Cycle 2 Cycle 3 Cycle 4 Cycle 5
15 1 9.96 8.77 11.79 27.46 21.13 2 3 3 3 2
16 1 14.72 44.09 24.53 12.47 15.57 2 3 2 2 2
19 1 0.00 0.00 0.16 4.40 6.88 0 1 3 4 3
20 1 0.00 0.00 0.00 0.00 0.13 0 0 1 1 1
63 1 8.08 4.95 23.62 8.02 5.63 2 2 2 4 4
64 1 2.25 1.34 18.18 5.35 5.54 2 2 2 4 4
72 1 6.29 0.99 15.72 44.37 20.24 3 4 4 4 4
73 1 0.00 0.01 0.00 0.02 0.01 0 1 0 1 1
129 2 0.97 33.67 9.19 12.44 - 2 2 3 4 -
131 2 10.76 25.21 22.73 12.43 - 3 3 3 3 -
133 2 0.00 0.01 0.00 0.63 - 0 1 1 3 -
134 2 0.00 0.00 0.00 0.00 - 0 1 0 0 -
148 2 2.58 12.14 6.52 19.33 - 4 4 5 4 -
153 2 0.00 0.00 0.17 0.03 - 0 1 2 1 -
154 2 0.00 0.00 0.00 0.00 - 0 0 1 0 -
155 2 0.19 0.04 0.05 0.01 - 2 4 3 2 -
203 3 0.00 0.00 0.00 0.00 - 0 0 0 0 -
207 3 0.15 0.90 3.82 0.46 - 1 2 2 1 -
215 3 3.81 10.11 4.42 0.03 - 6 5 5 2 -
219 3 0.14 0.02 0.02 - - 2 1 1 - -
226 3 0.00 0.31 0.01 0.01 - 0 1 1 1 -
230 3 0.16 0.32 0.27 0.04 - 2 2 2 2 -
233 3 0.00 0.02 0.12 0.29 - 1 1 1 2 -
238 3 0.30 1.40 13.58 5.30 - 5 3 4 3 -
303 4 0.00 0.00 0.00 0.00 - 0 0 0 0 -
304 4 0.00 0.00 0.00 0.00 - 0 0 0 0 -
306 4 3.41 4.03 1.62 0.29 - 3 4 4 2 -
307 4 0.03 0.04 0.01 4.75 - 1 2 2 4 -
309 4 33.48 55.72 65.03 38.95 - 4 3 3 6 -
311 4 0.00 0.00 0.00 0.00 - 0 0 0 0 -



Table 2 (cont.). Average % cover and number of invasive species detected in each plot by cycle. Cycle 1: 2007 – 2010; Cycle 2: 2011 – 2014; Cycle 3: 2015 – 2018; Cycle 4: 2019 – 2022; Cycle 5: 2023. Plots are sampled in 4-year rotations, with a quarter of the plots sampled each year (i.e. a panel). Note that cycle 5 is only partially complete. Plots with at least 10% average cover or 4 or more invasive species present are in yellow. Average cover and species counts may be increasing over time due to protocol changes, including additions to the indicator species list over time, a plot search for indicator species starting in 2009 and all woody species being sampled in quadrats starting in 2019. More details can be found in the MIDN Summary of Major Protocol Changes

table output
Average % Cover
Number of Species per Plot
PlotCode PanelCode Cycle 1 Cycle 2 Cycle 3 Cycle 4 Cycle 5 Cycle 1 Cycle 2 Cycle 3 Cycle 4 Cycle 5
314 4 0.01 0.02 0 0.00 - 1 1 1 1 -
317 4 0.00 0.00 0 0.01 - 0 1 1 1 -





Table 3. Number of plots out of 31 where priority invasive species were detected. Cycle 1: 2007 – 2010; Cycle 2: 2011 – 2014; Cycle 3: 2015 – 2018; Cycle 4: 2019 – 2022; Cycle 5: 2023. Species on the list include all exotic vascular plant species that are capable of dominating Eastern US forest habitats, and does not include all exotic species detected in MIDN forest plots. Note that cycle 5 is only partially complete.Average cover and species counts may be increasing over time due to additions to the protocols, including a plot search for additional indicator species starting in 2009 and all woody species being sampled in quadrats starting in 2019. Species have also been added to the indicator list over time. More details can be found in the MIDN Summary of Major Protocol Changes

table output
Latin Name Common Name Cycle 1 Cycle 2 Cycle 3 Cycle 4 Cycle 5
Ailanthus altissima Tree of heaven 3 1 0 2 0
Celastrus orbiculatus oriental bittersweet 1 1 1 3 0
Euonymus burningbush 0 0 0 0 1
Ligustrum privet 1 4 7 13 4
Lonicera - Exotic honeysuckle - exotic 3 1 0 0 0
Lonicera japonica Japanese honeysuckle 17 21 20 20 6
Microstegium vimineum Japanese stiltgrass 13 19 21 18 6
Persicaria longiseta Oriental lady's thumb 6 8 9 7 4
Persicaria perfoliata Asiatic tearthumb 0 1 1 0 0
Rosa multiflora multiflora rose 3 2 3 2 0



Table 4. Invasive plant and pest/pathogen early detections in Richmond National Battlefield Park observed in the most recent census of each plot. Coordinates are in UTM NAD83 Zone 18N.

table output
Plot Year X Y Latin Name Common Name Type
RICH-020 2023 297600 4163048 Litylenchus crenatae ssp. Mcannii complex beech leaf disease pest